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MEMBRANE ANALYSIS OF SCALLOPED SHELLS 


By Ravindra K. Vyas 
Ames Research Center 


SUMMARY 


The main objective of this report is to contribute to the understanding 
of membrane forces in scalloped shells, particularly scalloped paraboloids. 

It is shown how the stress function approach, used successfully in the mem- 
brane analysis of elliptic and hyperbolic paraboloids, can be applied to the 
scalloped paraboloids. In order to demonstrate the technique, two scalloped 
structures having practical significance as structural models for parachutes 
are discussed in detail. In the cases considered, it is shown that the stress 
function must satisfy a linear second-order partial-differential equation 
with variable coefficients. The basic method for solving this equation is 
given in the first appendix. For comparison purposes, numerical results for 
similar shells of revolution are also presented. These results show that the 
hoop and meridional stresses in the parachute type structure are lower than 
those in a corresponding shell of revolution. More significantly, the merid- 
ional stress in the parachute type shell decreases remarkably, and even 
becomes compressive, as the bulge of the scallop, defined by a positive param- 
eter 3, is increased. It is shown how a critical value, P cr , of this bulge 
parameter can be obtained below which the shell would remain in a compression- 
less state. A singular solution used in prescribing rope tension near the 
opening of the parachute type shell may also be useful in other aspects of 
parachute analysis. 

A brief review of the scope and the limitations of the method and the 
possibility of extending it to other scalloped shells of the same family, is 
presented in a separate section. Suggestions for further work are included 
in the concluding remarks. 


INTRODUCTION 


The advent of the space age is largely responsible for a growing inter- 
est in new structural concepts and techniques. Promising among the newer 
concepts are those of filament -wound structures and expandable structures. 
With these advances in technology, there exists a parallel need for theoreti- 
cal investigations of different structural forms. The membrane analysis of 
one such form, the scalloped shell, is presented in this report. 

As opposed to shells of revolution, scalloped shells have been rela- 
tively unexplored as structural forms . The parachute is , perhaps , the only 
form of scalloped structure that has been studied in any detail. An inter- 
esting history of the parachute and its development as a decelerating device 
is given in reference 1, which contains, in addition, valuable information 



about various aspects of parachute performance and design. The earliest pub- 
lished attempts at parachute stress analysis date from 1923 when the work of 
Taylor, Southwell, Griffiths, Jones, and Williams, on ideal shapes of para- 
chutes, was compiled by R. Jones in reference 2. Since then, several authors 
have continued this work, studying the many different aspects of parachute 
performance and design. 

Most of the work concerning parachutes has been based on certain simpli- 
fying assumptions regarding the character of stresses in the surface. These 
assumptions can be broadly classified into three categories as follows: 

(a) The circumferential stress in the parachute is either negligible (ref. 2) 
or negligible within a certain critical radius (ref. 3)* This position was 
taken by the early investigators; (b) It is the meridional stress which is 
negligible rather than the hoop stress (refs, k and 5)j (c) The circumferen- 

tial stress is proportional to the product of the applied pressure and a ref- 
erence radius. The meridional stress is assumed either negligible, for 
certain cases, or obtainable from the conditions of equilibrium (ref. 6). 

The present work differs from parachute stress analysis in two important 
respects. First, the emphasis here is on the membrane analysis of scalloped 
surfaces as shell structures. The shells discussed in detail belong to the 
general class of surfaces designated as bulged or scalloped paraboloids, of 
which the parachute type structure is only one member. Secondly, since the 
investigation is based on the membrane theory of shells, no a priori assump- 
tions are necessary regarding the nature of stresses - and none are made. 

The method of analysis relies on the stress function approach postulated 
in the membrane theory of shells of arbitrary shape. This stress function 
approach has been successfully used in the past for solving the problems of 
elliptic and hyperbolic paraboloids, and it is shown in this report that this 
method can also be successfully applied to a more general class of scalloped 
shells, which includes the scalloped paraboloids as special cases. In order 
to demonstrate the method, the solutions for two different types of scalloped 
structures are discussed in detail. Since the loading and boundary condi- 
tions prescribed for these two structures could also apply to comparable 
parachutes, it is conceivable that these structures could serve as structural 
models for parachutes. 

The technique used in this paper can also be applied to a wider variety 
of loading and edge conditions than those considered in the two examples. A 
more detailed discussion of the scope and the limitations of the method is 
given in a separate section. 


NOTATION 


A surface parameter used to define the scalloped paraboloids 
A m undetermined coefficients associated with the homogeneous solutions THm 
a reference radius, also maximum radius 
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a m,n 

b 

F(p,x) 

F m (x) ,F a (x) ,1 

Fp(x),F 7 (x) J 

f m ( x ) 

f n 

fa( x )^ f p( x ) 

g m (x),g a (x)/l 
gp(x) J 

h 

k 

k cr 

L 

m 

M 

N r ^-^0 .?^r0 

N r , Nq , N r 0 
n 

n r > n 0 J n r e 

P 

P 

Pr >P0 > Pz 


coefficients cf a power series in x (appendix B) 

undetermined coefficients associated with the homogeneous 
solutions ¥ 3 ^ 

radial distance to the inner opening 

coefficients of power series in x (appendix A) 

f met ion of p and x; denotes right hand side of an 
inhomogeneous partial differential equation 

functions of x 

f met ion of x 

normal component of the forces acting on the rope (fig. 9 ) 
functions of x 

functions of x 

normal distance from the apex of the shell to a horizontal 
plane passing through the point (r = a, 0 = 0 ) 

constant associated with a singular solution; also a subscript 

critical value of the constant k 

linear differential operator 

summation index; also subscript 

number of segments or gores 

membrane forces per unit length in the tangential plane of 
the shell (fig. 2 ) 

projected membrane forces (fig. 2 ) 

summation index; also denotes a power of r and powers of x 
dimensionless membrane forces 

concentrated load at the apex in the Z direction 
reference load intensity, also uniform normal pressure 
r, 0 ,z components of the load intensity 
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Q t >Vq>V z reduced load intensities 

Q(x) function of x (appendix B) 


9 lz 

R(x) 


s 


dimensionless load intensities 
function of x (appendix B) 

polar distance in cylindrical coordinates (fig. l), also a 
subscript 

summation index 


T 

u m (x) 

X 


y(x) 

z(r,e) 


rope tension 

function of x (appendix A) 

transformed coordinate corresponding to the angular distance 0 
in polar coordinates 

constant used in defining the range of the variable x 
function of x (appendix B) 

vertical coordinate defined as a function of polar coordinates 
r ,9 to define the middle surface of the shell 


a angle between a tangent to the middle surface and the r direc- 

tion; also used as a subscript and to denote a power of the 
variable p 

i ” 

a h^ a n^ a n coefficients of powers of x in the series solutions 
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bulge parameter used in defining the scalloped paraboloids; also 
used as a subscript 


cr 




7 

7 n’ 7 n 

6 

Sn 

9 


critical value of the bulge parameter 

coefficients of powers of x in the series solutions 

angle between a tangent to the middle surface and the Q direc- 
tion; also used as a subscript ^ denotes power of p 

coefficients of powers of x (appendix A) 

denotes a power of p 

coefficients of powers of x (appendix B) 

angular distance in cylindrical coordinates (fig. l) 
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A k 

A m 

P 

Pn 

i 

p 

Po 

®(r,0) 

¥(p,x) 

% 


H 

Hm 

P 


constant used in defining the range of the variable Q 
used to denote powers of p in the separated solution 
eigenvalues 

constant, describes the ratio of height h to the maximum radius a 
coefficients of power series in x (appendix B) 
constants of integration (appendix A) 

constant related to prescribed edge forces (fig. 10(b)) 
dimensionless radial distance 

dimensionless radius corresponding to an opening at r = b 
postulated stress function 

homogeneous solution corresponding to the separated form p^q)^(x) 
dimensionless stress function 

solution of a homogeneous differential equation subject to 
homogeneous boundary conditions 

homogeneous solutions corresponding to eigenvalues A m 

particular solution of an inhomogeneous differential equation 
subject to inhomogeneous boundary conditions 

Subscripts 

critical values of parameters 
homogeneous solution 

individual components of homogeneous solution 
particular solution 

MEMBRANE EQUATIONS 


Consider a thin shell whose middle surface is defined by the relationship 

Z = Z(r,e) 
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which expresses the vertical coordinate Z of the middle surface as a func- 
tion of the polar coordinates, r ,0 (fig* l) • A differential element of this 
surface is shown enlarged in figure 2 . The element is contained between two 
adjacent vertical planes, 0 = constant, and two adjacent vertical cylinders, 
r = constant. It is acted upon by the loads p r , p0 , p z (per unit area of 
the middle surface) and supported by the membrane forces N r , Mg, N r 0 • The 
latter quantities are forces per unit length of the line element through 
which they are transmitted . 


The membrane equations, which N r , % , N r 0 must satisfy in order to 
assure equilibrium,_can__be conveniently written in terms of the projected 1 
stress resultants N r , J% , N r0 and the reduced load intensities p r , p Q , p z 
as follows : 


i ± (r5r , + i ± (Sre , - h . - fr 


(la) 


i — ( r — H 

r dr \ dr r 


where 


1 d ,r= 


+ i (|z 5 ) + i 
r Sr \Se V r 


w r = 




P,. = 


Pe 


r Se 


a / 

dz 

Se \ 

( d r 

cos 

a 

cos 

~7 

cos 

J_ 

cos 

a 

S3 

hi 

CD 


cos 2 

a 


+ 7- OW 

or 


+ 2 


N r 0 


= -Vc 


*r 0 


r Se \ r Se 0 


(lb) 


(lc) 


n t 




(2) 


1/2 




cos a cos 7 r 

(cos 2 7 + cos 2 a sin 2 y ) 1 ^ 2 


cos a cos 7 


Pe 


_ (cos 2 a + cos 2 a sin 2 y ) 1 ^ 2 / 

P* = I ■ — ■ ■ 1 ^ - V- ( 


( 3 ) 


cos a cos 7 


tan a = — 

Sr 

tan r - i ^5 

r Se 


•‘•Projected on a plane z = const, (fig. 2 ) . 
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I 


With the introduction of the following dimensionless stress resultants 
and load intensities. 


n~ = 


N r ^ 
pa 


% V 

n 0 pa ' 


n r0 = 


Nr0 

pa 

P d 
y r 

9-r d 


> 

P Z 

<lz = TT 




(5) 


where p is a reference load intensity and a is a reference radius s r, 
the equations (la) - (lc) can be written as follows: 


dr 


r be 


1 _d_ 

r de 


d_ 

dr 


"d a Z + 1 dz 
,dr r dr. 


— n 


+ £ 


d 2 5 


n r0 + i 

y dr d0 r 2 d0 2 


4 - — 
r 


n r - n.0 
r 

ii 

i 

(6a) 


a 

(6b) 

dz ^ n r 

+ 1 dZ ^ n r0 


dr dr 

r d0 dr 


dz ^ n r0 

+ _1_ = _ 

(6c) 

dr d0 

r 2 de d0 a 


Equations (6a) - (6c) are applicable to any shell surface of arbitrary 
shape , in cylindrical coordinates , so long as there exists a one-to-one cor- 
respondence between points of the middle surface of the shell and their pro- 
jections on a horizontal plane (Z = const.)* If the forces n r ^ n q, n r0 are 
defined in terms of an Airy type stress function $(r ,0) such that 
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II 



( 7 ) 


n r 


l <1? + JL 

r dr r 2 be 2 




dr 

U 




*0 



n = _i_ ^ i d 2 fl> 

r 2 d0 r dr bd 


> 


then, it can be verified that equations (6a) and (6b) are automatically 
satisfied and equation (6c) reduces to: 


1 dZ + _1_ dfz\ dfj> + dZ _ 2_ d 2 Z d 2 $ 
r dr r2 dd 2 J dr 2 V 3 bd r2 dr 60/ dr d0 


+ 




d 2 0 + /l d 2 Z N j 6$ + / 2_ d 2 Z 
d0 2 \ r dr 2 / dr \r 3 dr d0 


_ 2 _ dz\ dl 

r 4 d0/ d0 


iz + dZ^£ + ldZ^£ 

a Or a r ^ q a 


, r oS 2 or/ J a 


H- 


8 2 Z 

dr 2 




(8 ) 2 


The governing equation (8) is a second-order partial differential equa- 
tion^ and, depending on the choice of surface, may have either constant or 
variable coeff icients . The stress function approach has been successfully 
used in the past for solving problems of elliptic and hyperbolic paraboloids . 
For each of these cases equation (8) has constant coefficients and is reduc- 
ible to a standard form: the former to the Poisson's equation, the latter to 

the wave equation. In the case of the scalloped paraboloids, equation (8) 
has variable coefficients, as will be seen later. 


THE SCALLOPED PARABOLOID 


The middle surface of the scalloped or bulged paraboloids may be defined 
by the relationship 

Z = Ar 2 (l + 00 2 ) , — 0 O £ 0 £ 0 O , Q q <. n , p > 0 (9) 3 


2 Reference 7 gives a similar formulation in rectangular coordinates . 

3 The shell surface, it may be noted, is symmetric with respect to the 

plane 0=0. 
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where A , 0 are constants. The parameter 0 defines the bulge of the 
scalloped surface and may he called the bulge parameter. It should be noted 
that when 0 equals zero , equation ( 9 ) describes a paraboloid of revolution, 
hence, the designation "bulged" or "scalloped" paraboloid. 

Such a shell, like any other structure, can be loaded in many different 
ways. As a roof structure, it may be loaded by gravitational and wind loads. 
As a pressure vessel bulkhead, it may be subjected to a uniform or varying 
normal pressure. As a parachute type structure, a uniform normal pressure 
would constitute a reasonable loading. In the analysis to follow, a uniform 
normal pressure will be assumed. 

Assuming the surface segment defined by equation ( 9 ), then the load 
intensities, the definitions ( 3 ) and (4) , and the equilibrium condition (8) 
take the form 


where 


Pr = 


sm a cos 7 




(cos 2 a + cos 2 a sin 2 7) 1/2 




cos a sin 7 


(cos 2 a + cos 2 a sin 2 7) 1/2 




cos a cos 7 


(cos 2 a + cos 2 a sin 2 7) 1/2 


Pr 

Vq 

P z 

lr 


% 

% 


tan a p 
tan 7 p 


-P 

tan a 
tan 7 

-1 


tan a = 2 Ar(l + 00 2 ) 
tan 7 = 2Ar 00 


( 10 ) 
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and 


(1 + 0 


+ pe 2 ) — + (1 + pe 2 ) — — - 2 pe 

Sr 2 r2 


1 d 2 $ 
r dr de 


+ (1 + P© 2 ) - — + 2p© • — — 

r dr r d 6 


_ 1 + Ar 2 

2Aa a 


U 



pe 2 + 4p 2 e 4 


(ii) 


Next, introduce the transformation 

p = r/a 
x =^P 0 

where a is a reference radius, and define 

<J>(r,0) = a a Y(p,x) 

h = pa 
p = constant 

The coordinate p in equations (12) is a dimensionless radius; the coordinate 
x, which is linearly related to 0, is introduced for convenience in alge- 
braic and numerical work, and the stress function ¥ in equation ( 13 ) is a 
dimensionless function. The height h, designated as the normal distance 
from the apex to a horizontal plane passing through the point r = a ,0=0, 
is defined in equations (l4) as a certain multiple of the reference radius a. 
If equations (12) - (l4) are substituted into equation (ll), it can he shown 
that the stress function ¥(p,x) must satisfy the differential equation: 


( 12 ) 

( 13 ) 

(14) 


L(¥) = F(p,x) 

where L is the linear differential operator defined by 


L = I (1 + P + x 2 ) 2— + (1 + x 2 ) - 2. 
L dp 2 P c>p 

+ (1 + x 2 ) -% - 2Px ± 


d 2 


P 2 dx £ 


P dp dx 


+ 2)3xil 

P dxj 


and 


F(p,x) = — + pp : 


S +{3P +^)x 2 + 4x 


(15) 


( 16 ) 


(17) 


In terms of the dimensionless stress function I and the variables p,x the 
membrane forces n r , ng , n r e can now be expressed as follows: 
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n e 


do 2 




- ppx 




n. 


r0 


. /fl/i a? i a 2 y 

VP 2 dx P dp d: 


(18) 


Membrane analysis of the scalloped shell, therefore, reduces to the 
mathematical problem of solving equation (15) subject to appropriate boundary 
conditions. The equation, as can be seen, is a linear second-order partial 
differential equation with variable coefficients, and includes all the deriv- 
atives of the stress function ¥ up to the second order. It does not seem 
possible to reduce it to a known standard form, though it can be solved, as 
will be shown in appendix A, by the method of separation of variables. 

In the next two sections, solutions of equation ( 15 ) will be discussed 
in detail, along with numerical results, for two particular types of struc- 
tures, namely the single scallop and the parachute type structure. The lat- 
ter is so called because of its similarity to a parachute in surface geometry, 
loading and boundary conditions. In each case, the numerical results for a 
similar shell of revolution will also be given for comparison purposes. It 
will be seen that the basic mathematical approach outlined in appendix A 
applies equally to both cases and yields solutions that are very similar in 
form. 


THE SINGLE SCALLOP 


Consider the symmetric scalloped shell shown in figure 3. It will be 
assumed for this particular example that the edge members in the vertical 
planes 0 = ±0 q cannot resist any thrust normal to their planes, so that 
Nq = 0 along these edges. Furthermore, it will be assumed that N r = 0 along 
the edge r = a. It must be mentioned that with the technique to be pre- 
sented in this section, it is also possible to prescribe other distributions 
of the force N r along this boundary. 

Since the structure, as well as the loading, is symmetric about the 
vertical plane 0=0, it is sufficient to assume a solution symmetric in 0 
and confine the analysis to the domain 0 £ r ^ a, 0 ^ 0 s£ 0 O . In terms of 
the stress function ¥ and the variables p,x, the mathematical problem of 
solving equilibrium equation (15), subject to the boundary conditions 
N 0 (r,0 o ) = 0, N r ( a , 0 ) = 0, can be described as follows: 

L(f) = F(p,x) O^p^l; 0 £x£x o ( 19 a) 
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■Ill I II 


II II 


= VP 2 * 2 at x = x 0 

dp 


1 dr + _p dfy 

-P dp P 2 dx 2 



at p = 1 


( 19 b) 

( 19 c) 


The variables in equation (19a) can be separated by assuming solutions 

of the form p^ k cp^Cx), where the functions cp^(x) satisfy a linear ordinary 
second-order differential equation with variable coefficients and can be 
obtained in the form of an infinite series in even powers of x. By use of 
this approach, 4 the solution ¥(p,x) can be obtained in the form: 

¥(p,x) = ¥p(p,x) + %(p,x) (20) 

The function ¥p in (20) satisfies both the inhomogeneous differential 
equation ( 19 a) and the inhomogeneous boundary condition ( 19 b), and has the 
form 

¥p(p,x) = — p 2 o^x 11 + m-P 4 ^ Pn^ 11 (2l) 

n— 0,2,4^ • • • n— 0,2,4, • • • 


The homogeneous solution 


%( p,x) has the form 


%(p,x) 


I 


^■m %m ( P^) 


( 22 ) 


where each function %m> defined by 

%m = p Am F m (x) (23) 

satisfies both the homogeneous differential equation 

L(%m) = 0 (24) 

and the homogeneous boundary condition 

fs (W - 0 (25) 

dp 

The values A m are the eigenvalues of the system described by equations (24) 
and (25) • These eigenvalues occur, in the present case, as pairs of posi- 
tive and negative numbers of equal magnitude, so that the corresponding 

homogeneous solutions take on the form p m f m (x) , p g m (x) , where 
4 For more details see appendix A. 
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(26) 


f m (x) = ^ 


C m,n x 


n 


II 0;2;4^ • • • 


S m ( x ) = ^ %,n xI1 


n=o,2,4 , 


The homogeneous solution can now he written as 


% - 


A 




m 1 , 2 ; 3 ; • • • 


m 


n— ’O , 2 ) 4 , • • 


C m,n xn + 


-A 


Bo 

m 




■m 


n=o,2,4 , . 


n 


D m,n x 


( 27 ) 


Substituting equations (21) and ( 27 ) into equation (20), the solution 1 may 
he written 


¥(p,x) 



a n x 


.n 



n=o,2,4 


Pn x 


,n 


A 


A m P 


m= l , 2 , 3 , 


m 


Ii~Oy 2 y^ y . . . 


■ / m.n 


n 

x 




I 



D m,n x 


n 


( 28 ) 


m=i,2,3., . . . n=o ; 2,4 ; • . . 

Applying ( 28 ) and definitions (l 8 ) , the stress resultants can he expressed as 
follows : 

n r = ~ [ 2 a n + P(n + 2 )(n + l)a n 4 ^]x n 

n = 0,2 ,4 y . . . 

+ PP 2 ^ t^Pn + P( n + ^)(n + 1 )Pn+2^ x3:i " pp 2 (^— + x ^) 

n = 0 y2 ,4 , • • • 


+ ^ A m P ^ + P( n + 2) (n + l)C m ^ n+2 ]x n 

m= 1,2 y3 y . . . n=o , 2 , 4 ,... 

+ ^ m ^ t-A m D m ,n + P( n+ 2) (n + l)l ) mj , n +2 J xI1 (29) 

m=i,2,3 , . . . n-o,2,4, . . . 
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EH— 1^2 ? 3 , ••• H 0^2 ;4; ■ • • 


+ Y B mP (Am+2) X \n(\n + D\/ (30) 

m=i,2,3, . . . n=o,2,4,... 


n r0 = VP 


2ji 


I 


n-i 2 
ncx^x + (j.p 


I 


- 3nP n x 


n-i 


11 = 0,2 ,4 




n=0,2 ;4 V . 


I 


A mP 


Am-2 


n 


(l - A m )C m n x 


n-i 


HI — 1 ,2. }3 ) • • • El O ,2 ^4 y • • • 

-(A m +2) 


+ > B m p 

m= 1^2 ^3 ^ • 


n(l + A m )D m ^ n x 


(31) 


n=o .2 >4 


It should be observed that the homogeneous solutions associated with the 
negative eigenvalues are singular at the point p = 0, which is an included 
point of the do ma in ■ Consequently, the undetermined coefficients B m asso- 
ciated with these solutions must be chosen as zero. The coefficients A m 
can be chosen appropriately to prescribe suitable edge forces at p = 1 
(correspond ing to r = a) . In particular it is possible to choose these 
coefficients so as to satisfy the condition prescribed in (19c) • For parab- 
oloids of revolution (0 = 0) 5 the eigenfunctions f m (x) form a set of 
orthogonal trigonometric functions, and the coefficients A m can be con- 
veniently determined by means of this orthogonality property. In the present 
case (0 f 0), no such orthogonality property is readily available, and it is 
necessary to resort to suitable numerical methods such as collocation, match- 
ing coefficients of power series, or mean square approximation. The last 
technique has been found to give the most satisfactory results, and its 
application to the case of scalloped paraboloids is described briefly in 
appendix B. 


transformation (11 ) for this case is 9 = x. 
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Consider, as a particular example, the case 3 = 1, [± = 0.5, 0 o = rt/ 6 - 
The homogeneous solution (22) was found to be sufficiently approximated by 
the six-term truncated series 

6 6 

% = A m%m = Y, AmP m Z (32) 

m= 1,2 ; • • • m= 1,2 , • • • n=o ,2,4, • • • 

where the eigenvalues A m and the coefficients A m are listed in table I. 
These values have been obtained according to the methods discussed in the 
appendixes A and B. The coefficients o^, 3 n associated with the particular 
solution ( 2 l) and the coefficients C m ^ n associated with the homogeneous 
solutions i farn may then be evaluated by means of the recursion relationships 
given in appendix A. In this example, the first 25 terms of each of the 
series expansions in x in equations ( 29 ) - ( 31 ) were used to obtain the 
stress resultants presented in figures 4-7. 

TABLE I.- SINGLE SCALLOP: EIGENVALUES |A m | AND THE 

UNDETERMINED COEFFICIENTS A m 


m 

l*ml 

A m 

1 

2 .26514 

0.453853 

2 

6.63706 

-.913284xio -2 

3 

11.0383 

.689597X10 -3 

4 

15-4446 

- . 121309 x 10 _3 

5 

19.8525 

. 134767 x 10-4 

6 

24.2612 

- . 286 o 88 xio _4 


The broken line in figure 4 represents the meridional force n r at the 
edge p = 1 due to the particular solution ¥p; the solid line represents the 
force n r due to this particular solution plus the homogeneous solution 
( ec l* (32)). This homogeneous solution is added to obtain the desired 
vanishing edge force n r along the edge p = 1. The difference between the 
broken and solid lines in figure 4 indicates the efficiency of this edge 
correction by the six-term approximation (32) . This efficiency would improve 
as more solutions were included in (32). The forces n© , n r e along 

this edge are plotted in figure 5 - 

For the interior points of the domain, it can be shown that only the 
first few homogeneous solutions are of importance, since the solutions cor- 
responding to higher eigenvalues rapidly diminish as p decreases from 
p = 1. The solid lines in figures 6 and 7 represent the stress resultants 
in the scallop, while the broken lines represent the stress resultants in the 
corresponding shell of revolution (3 = 0) . The forces n r , n@ along the 
central meridian (fig. 6 ) are lower in most parts of the scalloped shell than 
in the shell of revolution. The same is true for shear stress along the 
edge Q = Q Q9 although along this edge the meridional force n r is higher in 
the scalloped shell for a considerable distance from the apex. 
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It should be observed that for both types of shells the stresses tend to 
be relatively high in the vicinity of the three corners. These stress con- 
centrations are pronounced, in the case of the scalloped shells, near the 
apex, and in the case of the shells of revolution, near the lower two corners. 
For either type of structure the corner regions need special attention and 
should be designed to withstand higher stresses than the rest of the shell. 

The method of analysis for a shell terminated at an inner boundary 
r=b/0(p=p o ^0) would be similar. In this event, however, the solu- 
tions corresponding to the negative eigenvalues would have to be included in 
the analysis, since these solutions are necessary for prescribing appropriate 
boundary values at r = b. If the lowest eigenvalue |A 1 | is sufficiently 
high, then the coefficients Am; B m associated with the two sets of homoge- 
neous solutions can be determined separately, by mean square approximation at 
each of the two boundaries. Otherwise, it may be necessary to use an itera- 
tive scheme, similar to the one described in reference 8 , in order to obtain 
a satisfactory solution. 

It should be noted at this point that the bulge parameter P has been 
prescribed unity only for the purpose of the numerical example . The tech- 
nique places no restriction on the magnitude of this parameter, although, of 
course, it must always be positive, since surfaces with negative P belong 
to a different family. 


PARACHUTE TYPE STRUCTURE 


The method demonstrated in the preceding section can also be applied to 
the parachute type structure with an opening at the top, as shown in fig- 
ure 8 . In this case, however, the circumferential force n^ need no longer 
be zero at 0 = ±0 Q . Because of structural symmetry and symmetric loading it 
is sufficient to solve the problem, as for the single scallop, in the domain 
b ^ r £ a, 0^0 ^ 0 q • 


The boundary condition at 0 = 0 Q can be obtained as follows. If the 
rope is considered an individual structure, loaded by the membrane forces 
from the shell at the edge 0 = 0 O , it can be shown that the rope tension T, 
at any coordinate r, must satisfy the equilibrium condition (see fig. 9 ), 


T = 


[1 + 4 A 2 r 2 ( 1 + P 0 2 ) 2 ] 372 

. f n 

2A(1 + P© 2 ) 


4Apap0 o r[l + lJA 2 r 2 (l + (30 2 ) 2 ] 1/2 
2A(1 + pej) 


( 33 ) 


On the other hand, considering the entire structure above a given radius r , 
it follows that the rope tension T and the membrane forces n r , n r g must 
also satisfy the condition of equilibrium of vertical forces: 
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p[l + 4A 2 r 2 (l + (30 2 ) 2 ] 

jtr 2 

P_ 

2A(1 + f30 2 )r 

N 

" pN 


r°o 

p®o 

- 4aAr 2 / (l + (B0 2 )n r d0 - 4Aar 2 

/ 00n 0 d0 

J o v. 

' 0 


(3*0 


where P is an arbitrary vertical load concentrated downward at the apex, 
and N is the number of scallops. The load P may be assigned any value 
appropriate to the problem of interest. Eliminating the rope tension T 
between equations ( 33 ) and ( 3 * 0 , the boundary condition at 0 = 0 O is 
obtained: 


r 2 (30 o n 0 (r,0 o ) 



+ pe 2 )n r 


+ (30n r0 ]ae 


jit 2 _ P 
4AEa 4ApaN 


(35) 


The boundary conditions at the edges r = b and r = a depend on the 
nature of the structure . For a model parachute it would be reasonable , for 
example, to expect an almost vanishing meridional force n r along these 
boundaries: such a surface would be expected to carry the applied loads by 

means of tensile and shear forces alone, since the material may be assumed 
incapable of sustaining compressive stresses. Assuming uniform normal pres- 
sure, the governing differential equation, condition ( 35 ) and the above 
requirements can be described as follows: 


L(¥) = F(p,x) p Q £ p s: 1 ; 0 <; x £ x Q 


(36) 



x=x Q 





JL 

p 2 dx 2 / 


+ 


P: 


_1_ dT 

P 2 dx 


1 d a ¥ \ 
P dp d yj_ 


dx 


x oP 2 

r 



Pn/P 

4pa 2 pN 


(37) 


n r « 0 
n r « 0 
n r , 

The solution of equations 
can be written in the form: 


at p = p 0 (r = b) 

(38) 

at p = 1 (r = a) 

(39) 

n 0 s 0 

{bo) 


( 36 ), by the method described in appendix A, 
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¥(p,x) = i p 2 a n x n - A log(p/l + x 2 ) 

n=0,2,4, . • . 

+ PP 4 £ P n xI1 + ^ A mp Am Yi Cm ’ n ^ 

n=0,2,4 } ... m=l,2,3,... 11=0,2,4,... 

+ Y EmP‘ Am X D m,n^ 

HI - * 1 j2 j3 J • • • H O y 2 } 4 j • • • 

= *p + % 

where 

= g Y a n^ + ^ ^ ^ ^ lo s(P h + 

n-O ^2^4^ • • • n o ^ 2 ^4 ^ • • • 

B mP Am Y 

m- 1,?2^3^ • * • n— o j2j4j • • • m — i ^ 2 ^ 3 ^ • • • n = o y 2 y 4 y 



k . ^ r i 1 

2pa 2 N L(1 - Xq)P + p tan 1 x Q + (xq/3)_ 


(4l) 


m,n 


XH 


(42) 


A comparison of equations (4l) with (28) rev eals th at the forms of the 
two solutions differ only in the term k/2 log ( p \J 1 + x 2 ) . By means of 
equations (4l) , the membrane forces n r , n q, n r @ and the rope tension T can 
be expressed as follows: 


n~ = - 


_k_ _1_ 

2p p 2 


P(1 - x 2 )- 1 


1 + 


2p 


[20^ + P(n + 2)(n + l)o- n+2 ]x n 


n=o,2,4 ; ... 


+ [iff 


n 0 < y2^4^»»» 


(1 + X 2 ) 2 J 

Y ^4p n + 3(n + 2)(n + l)3 n+2 ]x rl “ (j + x2 ) 

A mP Am_2 Y [AmC ^ n + P(n + 2)(n + l)C m,n +2 J xn 

?2 y3 y • • • n-O, 2 ,4 , • • • 

^ t B m p m ^ L“^m^m,n + ^( n + 2)(n + l)^m^n-l- 2 -l x ' 


HI— 1 , 2 j 3 j • • • 


n-O ;2 ;4 ; • • • 


(43) 


l8 


I 



I 


ne 


_k_ _L 

2n P 2 


2 n 


I 

n=0,2,4. 






* I 


An/”' 2 I 

11=0,2,4, 


\i(\i ~ -*-) < - : m,n xrL 


+ £ B m p' (Am+a) £ A m (A m + l)D m ^ n x n (44) 


n r 0 \/"^ 


_k_ JL_ x 
2 p p 2 1 + x 2 



n=o ; 2 ; 4 


n-i 

na n x 


- M-P 2 


I 


11 = 0 , 2,4 




Am -2 


+ ) A m p 

in — 1,2,3, • • • n o ,2,4, . * • 


+ ) B m p 

HI— 1,2,3 , • • • 


^ n(l - A m )C m ^ ri x n 

-(Am«) £ n(1 + Am)Dm<nX n 


n— O ,2,4, • • • 


(45) 


T = 2 pa 2 \/0 x Q p 


[1 + 4|i2p2(l + xg) 2 ] 1/2 

(1 + *g) 


ne(p^x 0 ) 


(46) 


The boundary conditions ( 38 ) and (39) and the requirement (4o) represent 
perhaps an ideal situation, but the tensile requirement (4o) could be satis- 
fied even if the forces n r were tensile at the upper and lower edges. 


The stresses due to the particular solution are of interest in view 

of requirement (4o) . In the present investigation it has been found that for 
a given choice of the parameters N, p Q («l) and p the above requirement 
will be satisfactorily met if (3 and k are less than certain readily deter- 
mined critical values |B cr and k cr . The significance of these parameters 

will be discussed later. It should be observed in passing that for practical 
parachute models N is generally large (say >10). In the present study it 
has been found that values of |A m | depend largely on N and are of the order 
of magnitude mN. For large N, the particular solution then becomes 

the dominant solution, and the approximate solution ¥ « can be used to 

obtain satisfactory values for stresses in most parts of the shell. Consider, 

as a numerical example, a parachute model with 12 scallops (N = 12 ), p Q = 0 .1 
and p = 0.4. For this combination of the geometric parameters, the critical 
values of (3, k have been determined to be 0.5 and 0.0017,? respectively, and 
the corresponding solution ¥p meets requirements ( 36 ) - ( 38 ) . It may be 
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recalled here that, as for the single scallop, the two sets of homogeneous 
solutions associated with the positive and negative eigenvalues are to he 
used to satisfy proper edge conditions at the lower and upper boundaries. 

Since condition ( 38 ) at the upper edge is already satisfactorily met, the 

coefficients B m , associated with the solutions in p ^ m , may he appropriately 
chosen as zero. Furthermore, with P = P cr , the edge forces n r at p = 1 
due to the particular solution are already very small, and a five-term trun- 
cated series for % satisfactorily meets requirement ( 39 ) . The correspond- 
ing five values of A m , A m are listed in table II. 

TABLE II.- PARACHUTE TYPE STRUCTURE: THE CONSTANT k, 

THE EIGENVALUES |A m | AND THE UNDETERMINED 
COEFFICIENTS A m 

k = 0.001726 


m 


Am. 

1 

12.7379 

- 0 .248942X10 -3 

2 

23-3730 

. 200055 XIO - 4 

3 

33-5535 

- . 488396 XLO -5 

k 

43.5804 

. 179266 KLO - 5 

5 

53.5374 

- .115634X10 - 5 


Before discussing the results for this case and a similar shell of revo- 
lution, it is necessary to contrast certain features of the solutions for the 
two types of structures . For both the scalloped shell and the shell of revo- 
lution it is possible to find a solution ¥ p that meets requirements ( 36 )- 
( 38 ) . In the shell of revolution, the functions fm form a set of orthogonal 
trigonometric functions. The coefficients A m can then be readily deter- 
mined, from this orthogonality property, to prescribe appropriate edge forces 
at p = 1. Along this edge the force n r due to the particular solution Ip 
is equal to (l - and the dimensionless sum of its vertical components 

along this edge is a downward force of magnitude 0 o (l - Pq) • The sum of the 
vertical components of the force n r due to the homogeneous solution 
vanishes along this edge. This implies that a proper boundary condition at 
p = 1 consists of a distribution of the force n r , which has the same verti- 
cal resultant as that due to the particular solution. Two such distributions 
are shown in figures 10 (a) and 10(b). For the scalloped model, on the other 
hand, it is not necessary to follow the same pattern because the correspond- 
ing homogeneous solutions do contribute a vertical resultant. In this case, 
however, the coefficients A m have to be determined by numerical techniques, 
with resulting numerical difficulties if N is large and if a large number 
of terms of the homogeneous solution are needed to prescribe vanishing edge 
force n r at p = 1. However, if P is close to Per; then only a few terms 
of the homogeneous solution would be adequate to meet requirement (39)* This 
is true because in the scalloped model, as will be seen later, the forces n r 
decrease as p increases and become very small at the base when P is close 
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to 0 cr . In the shell of revolution, on the other hand, these forces increase 
toward the base and become almost constant in the lower parts of the shell. 

The membrane forces n r , ng, n r Q and the rope tension T presented in 
figures Il-l4 have been obtained by means of the first 25 terms of each of the 
series in x in expressions (43) - (46) . The broken lines in these figures 
represent the membrane forces in a similar shell of revolution subjected to 
the boundary condition in figure 10(b) for | = 0.1. It can be seen from 
these figures that, for equilibrium of the shell, all three stress resultants 
n r , n Q > n r0 are necessary. For the particular parameters chosen the force 
n r decreases toward the edge p = 1. This drop in n r is more noticeable 
along the central meridian than along the edge 0 = 0 O (fig. 12). The hoop 
force tiq at a given radius (fig. 13 ) remains nearly constant although its 
magnitude varies with the radius and reaches a peak value at the edge p = 1 . 
The shear forces are relatively low, except in the lower half of the shell 
near the edge 0 = 0 Q . The rope tension T (fig. l4) has a maximum value at 
the base (p = l) and decreases steadily toward the apex. 

It is interesting to note how the membrane forces in the scalloped shell 
differ from those in a similar shell of revolution (figs. 11 and 12). It can 
be seen that the forces n r in the scalloped shell are significantly lower 
than those in the shell of revolution. The forces n@ are also lower but not 
as significantly as the forces n r - The shear force n t q in the shell of 
revolution vanishes at 0 = 0 Q , whereas in the scalloped shell a relatively 
low but nonvanishing shear force is necessary at this edge. The behavior of 
the forces n r , n^ in the shell of revolution along 0 = 0 Q is also worth 
noting: These forces have high magnitudes in the neighborhood of the lower 

corners, and ng even changes sign and becomes compressive. 

It sh ould be mentioned, at this point, that the solution k/2p log 
(Psjl + x 2 ), with k = 0.0017 in the particular example, corresponds to a down- 
ward load of magnitude 0.35tfb 2 p at the apex. In figure 8 , the ropes, but 
not the scallops, are continued to the apex. Prescribing a value for k 
therefore amounts to assuming a certain amount of tension in these ropes at 
p = p Q . In the particular example considered, if k > 0.0017 the edge forces 
n r at p = p Q due to Tp would be compressive, while if k < 0.0017 these 
forces would be tensile. In either case, appropriate values for n r would 
have to be prescribed along this edge by means of the homogeneous solutions 
associated with the negative eigenvalues. 

It was mentioned earlier in this section that in order to meet the 
tensile requirement (4o), the parameter 0 should be less than 0 cr * This 
can be illustrated in figure 15 where membrane forces due to the dominant 
solution ¥p are plotted along the central meridian for different values of 
0. The constant k in each case has been chosen equal to k cr to meet 
requirement ( 38 ) . This figure shows that when 0 is very small the membrane 
forces approach those in a corresponding shell of revolution (0 = 0) . How- 
ever, as 0 increases the stress resultant n r decreases markedly, and, for 
0 = 0 . 6 , even changes sign and becomes compressive in the region p ^ 0 . 75 * 

The hoop force n q decreases also, although not as significantly as the 
meridional force n r . It can be seen that when 0 is less than 0.5j the 
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shell (and a corresponding parachute) could support the applied uniform pres- 
sure by means of tensile and shear forces alone. 


SCOPE AND LIMITATIONS 


The stress function approach is valid as long as there exists a one-to- 
one correspondence between points of the middle surface and their projections 
on a reference plane. The scalloped paraboloids (eg.. ( 9 )) meet this 
requirement. The mathematical approach in appendix A for finding this stress 
function f relies on two important conditions. First, the differential 
equation and the boundary conditions must be such that the method of separa- 
tion of variables can be applied. Second, the solutions involving the vari- 
able x must be obtainable in the form of an infinite series in powers of x. 
These requirements are fulfilled in the problems considered so far. It is of 
interest to mention other problems where these requirements may be met . 

Consider, for example, the more general class of surfaces 

Z = Ar n (l + 00 2 ) (47) 

of which the scalloped paraboloids constitute a special case. These surfaces 
meet the above requirements for application of the stress function approach. 
When such shells are subjected to uniform normal pressure, the stress function 
¥ must satisfy the differential equation. 


U + 0 + S J] Sfl _ 2 (n - l)(3x i -ifl- + i n(n - l)(l 
V 2 2 J ap 2 P dp dx 2 


+ x2) 4 

P 2 dx 2 


+ - n(n - 1) (1 + x 2 ) i + 2(n - l)/3x -i- 


p a P 


2 -n 


p 2 dx 


1 ^ “ , n 

= — P + PP 

2p. 


— — — + fn 2 — — - + Sp'jx 2 + n 2 x 4 
jci + 1 \ n + 1 / _ 


(48) 


where the membrane forces are related to the stress function 1 as follows: 


i jl d 2 1 

> c 

d 2 ¥ 


a r = 7f + 2 = " w 

p dp p 2 a x s 


n 


n + 1 


+ x 2 


o I n? 

„ _ rnfl dy 1 d 2 Y 

n r0 - n/P H r T- - - 


D 2 dx P dp dx. 


( 49 ) 
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It can "be verified that the boundary conditions , similar to those for the two 
cases considered, are such that the method of appendix A can be applied to 
equation (48) to obtain solution ¥ in the form: 


¥(p,x) = i 
2|JL 


a. 


'Fa(x) 


+ pp 7 Fp(x) + A 

2p 


P B F 7 (x) 


+ %(p,x) 


(50) 


where 


a = 4 - 


n 1 


7 = n + 2 ! 
& = 2 - n 


k 


= 0 for a single scallop 
£ 0 for a parachute model 


(51) 


Whe n 5 = 0 (or n = 2) the third term of (50) is replaced by k/2p log 
(Pn/ 1 + x 2 ) . 


It should be pointed out that it is not generally possible to find a 
solution in the form of equation ( 50 ) if any of the indices a, 7 , 6 is 
equal to an eigenvalue of the corresponding eigenvalue problem. For example, 
A = 1 happens to be an eigenvalue for all n; hence, difficulties may be 
encountered when n = 1 and n = 3* In order to predict whether a solution in 
the form of equation ( 50 ) does or does not exist for a given problem, it is 
necessary to investigate each individual case f urther . It is interesting to 
note that since the shells corresponding to n > 2 become relatively flat in 
the polar regions, the membrane forces must be high in these regions in order 
to balance the applied normal pressure. For such cases the structure would 
require additional reinforcement in the upper regions. 


CONCLUDING REMARKS 


In the preceding sections it has been shown that the stress function 
approach can be successfully used in solving problems of scalloped parabo- 
loids . It has been further demonstrated that the surface parameters associated 
with these problems may profoundly affect the application of the scalloped 
shell - particularly where a compressionless state is desired. The two 
scalloped structures considered in detail in this report are characterized by 
the parameters N, p, p Q , P, and defined by equation ( 9 ). It has been shown 
for the parachute model that for a given large N and given values of the 
parameters p , p, it is possible to use k = k cr in the particular solution 
to determine the maximum allowable bulge parameter P cr , "beyond which the 
meridional stress could become compressive. It was found in the particular 
example considered that for given N (=12) and p Q (= 0.1 ), the greater the 
value of p, the smaller the value of Per* lb would be useful to carry this 
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investigation further and plot relationships between p and P cr for differ- 
ent sets of p Q and N. Such plots could serve as a useful aid in the design 
of parachute -type structures . 

As mentioned earlier in this report, the only difference between the 
for ms of th e solutions for the two cases considered is the term k/ 2 p log 
( p + x 2 ) , which is singular at the apex. For the single scallop, this 
term is absent. For the parachute model shown in figure 8 this term is used 
to prescribe a certain amount of tension in the ropes in the region p < p Q . 
Solutions of this type may be useful in other aspects of parachute analysis, 
as well. 

In the case of the single scallop the character of the stresses does not 
seem to be radically affected by the magnitude of the bulge parameter P; but 
such is not the case with the parachute type structure, as can be seen In 
figure 15 * The latter structure, when designed to the critical p, has a 
lower hoop stress and a much lower meridional stress than the corresponding 
shell of revolution. In comparison with shells of revolution, indications 
are that the scalloped shell may be an efficient structural form warranting 
further study. It would be useful indeed to extend the analysis to include 
bending, deformation, and material properties, and to consider the scalloped 
shell for other structural applications such as pressure vessel bulkheads. 


Ames Research Center 
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APPENDIX A 


THE SOLUTION ¥ FOR THE SINGLE SCALLOP 
PARTICULAR SOLUTION 


Consider the differential equation for the single scallop, 
L[Uf(p,x)] = F(p,x) ; 0 s; p s; 1 , 0 <s x ss: x Q 

subject to the boundary condition 

(^| - MP 5 ^ 2 ^ = 0 at x = x Q 

The operator L and the function F are defined in equations (l6) anri 
Let 

¥(p,x) = ¥ p (p,x) + %(p,x) 

where ¥ p satisfies equation (Al) and the boundary condition (A2) . If 
is assumed in the separated form, 

¥p( p >x) = JL p 2 F a (x) + pp 4 Fp (x) 

then F a , F p must satisfy the equations 

P(l + x 2 )F^ - 2PxF^, + 2(2 + P + 2x 2 ) F a = 1 

P(l + x 2 )Fp - 6pxF p + 4(4 + 3P + 4x 2 )F p = | + ^3P + x 2 + 4x 4 

0 £ x s» x Q 


and the boundary conditions 


F a (x 0 ) = 0 
x o 

V*o) = 


x When used as function superscripts, the primes denote appropriate 
ferentiation with respect to x. 


(Al) 

(A2) 

(17). 

(A3) 

w 

P 

(a4) 

(A5) 1 

(A6) 

(A7) 

(A8) 

dif- 
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Solutions for F a (x ) , Fp(x) can be found in the form of a convergent 
series containing terms in even powers of x as follows: 

F a( x ) = f a( x ) + ^ a Sa( x ) 

Fp(x) = fp(x) + Ppgp(x) 
where p a , |ip are constants to he determined and 


f aW = Y, 


a n x 


.n 


0 , 2 , 4 , . 


fp(x) = y 


0 , 2,4 , . • . 


g a (x) = y ^ 


o ,2 , • 


gp(x) - 




n 


(A9) 

(A10) 

(All) 

(A12) 

(A13) 

(A14) 


0,2 ,4 


The functions f a , fp must satisfy the inhomogeneous equations (A5) and 
(A6) , while g a , gp must satisfy the corresponding homogeneous equations. 

The coefficients <x^, cx^, p^, P^ can then he obtained from the resulting 
recurrence relationships. 


(A15) 
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1 
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It 1 fl 
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for n = 0 , 2 , b , . . . 
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where 


and 
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The constants of integration M-p in (A9) and (A10) are obtained 

requiring that Fp satisfy conditions (A7) and (a8), so that 


Ha 


O , 2 . ,4 , ■ » . 




L 

0 , 2,4 , ■ 


ii n 

°h X g 


(xg/12) - E P^xg 
HP 

ft -p 

2 Pn x o 


Defining the coefficients a n , P n by 


“h = a n 

Pn = PA 


+ H a a n 



n = 0 , 2 , 4, . . . 


The solutions F a and Fp can be finally written in the form: 


F a (x) = £ a n x n 

n=o ^2 ^4 ) • • • 

Fp(x) = ^ C^x 11 

n=0 ,2 ,4 } . • . 


HOMOGENEOUS SOLUTION 


Consider a homogeneous solution 

= P Am f m(x) 


(A19) 

Dy 

(A20) 

(A2l) 

(A22) 

(A23) 

(A24) 

(A25) 
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which satisfies 


L (%J = 0 (A26) 

and the boundary conditions 

%m(p> x ) = 0 at x = x Q (A27) 

dp 

The function f m (x) , then, must satisfy the equation 

|3(1 + x2)f^ - 2f3(A m - l)xfi + A m [P(A m - l) + A m (l + x2)]f m = 0 (A28) 

subject to the condition 

A m (A m - l)f m (x c ) =0 at x = x Q (A29) 

A nontrivial solution for f m (x) exists only when A m equals an eigen- 
value of the problem defined by equations (A 28 ) and (A29) • The values A m = 0 
and A m = 1 are obviously two eigenvalues. However, the corresponding 
homogeneous solutions make no contribution to the membrane forces n r , n^, 
n r @ . The eigenvalues of interest, therefore, are associated with the 
condition 


f m U 0 ) = 0 (A30) 

It should be observed that the solution for f m , when one exists, can be 
obtained through the method of infinite series, that is, f m (x) may be expanded 
in the form, 


fmW 
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(A31) 


where 


c m,o = 1 


"m,n 


= 0 


for n < 0 


and 
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^m,n+2 


-1 

(n+l) (n+2) 
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n(n+l) - 2n(A m -l)+ fl + i 

1 ^m ~ \n 

Am 

C m,n + p c m,n-2 


- \ H/ 

— 



for n = 0 , 2 , 4, . . . 


(A32) 


Equation (A3l) can be applied to (A 30 ) to obtain an infinite series in powers 
of A m . The zeroes of this infinite series are the required eigenvalues. 
From a practical standpoint, the leading eigenvalues can be found by truncat- 
ing such a series after a finite number of terms and then finding the roots 
of the resulting polynomial in A m . In the present work, it was found that 
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for large N(= n/9 Q ) the polynomial approach gives convergent values only for 
the first one or two zeroes of the infinite series . A more satisfactory 
approach in this case has been found to he the method of interpolation or the 
false position method (ref. 9 ). It must he mentioned, however, that this 
method is to he applied to equation (A30) in a special way and not to the 
truncated polynomial. The process essentially consists of selecting a set of 
values for A m at chosen fixed intervals, obtaining the coefficients C m n 
from (A 32 ), and using these C m ^ n in (A3l) to evaluate each f m (x 0 ) . The 
approximate zeroes of (A29) are then easily located and it is only necessary 
to repeat the entire process a suitable number of times in order to converge 
to suitably accurate roots. In order to use this approach most effectively, 
it is advantageous to have prior information about the character of the eigen- 
values A m - In this connection the following discussion will be of interest. 

If, instead of using the form (A3l), f m (x) is expressed in the following 

form: 

f m (x) = (1 + x 2 ) 2 u m (x) (A33) 


then it can be verified that u m must satisfy 




+ 



(1 + x 2 ) 2 


A™ 

/v m 


%i = 0 


(A3*0 


%i(x 0 ) = 0 


(A35) 


Furthermore, if A m is a complex eigenvalue and u m the corresponding com- 
plex eigenfunction, then their complex conjugates A m , u m also satisfy 
equations (A3*0 and (A35) • 

Using u m , u m in Green’s first formula (ref. 10) 


f ° %% dx = (^1%)^° - / KlKi 


dx 


(A36) 


and realizing that u m (x 0 ) = 0 and (because of symmetry) u^(o) = 0 , equa- 
tion (A 36 ) reduces to 




KK dx 


Substituting for 




from equation (A 3*0 


and rearranging yields 


(A37) 
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, r» x o 

%% dx + J Q 


(l + X 2 ) 2 


u m u m dx 



1 

_(l + x 2 ) 2 


i 

p 


"nfta 


(A 38 ) 


where JB > 0 by definition . 


It is evident from the above expression that, since ^ is a nontrivial 
solution, A^ must always be positive definite. Hence, the eigenvalues are 
always real , occurring in pairs of positive and negative values of eq. ua! mag- 
nitudes . Consequently, the homogeneous solution % can be written as 


%(p,x) = 


I 


m-i, 2 ; 3 ; 


(A39) 


where the functions f m , g m correspond to the eigenvalues -f-A m and -A m , 
respectively. The function fm is obtainable in the form defined in 
equation (A31), and g m is obtainable in the form 


SmW 


z 

n=o, 2,4 


A21 ,n x 


(a 4 o) 


where the coefficients D m n can be generated from the recurrence relation 
(A32) by means of the negative value -A m * The stress function can then be 
written as 


¥(p,x) = -±- 

2[i 


ot^x 11 + pp 4 ^ 


(3 x J 
K n 


.n 


n=o, 2 ,4 


n O ,2 ,4 j • 


A mP' 


A: 


■m 


c m,n x 


n 


m— 1 y 2. j • • • n— • • * 


+ ^ B mp" Am X ( A ^) 

xn ~ i ) 2 ) 3 ) . , . 11=0,2,4, . . . 

The homogeneous solutions associated with the negative eigenvalues tend to 
infinity as p approaches zero. Hence, the coefficients associated 

with these solutions must be chosen to be zero if p = 0 is to be an included 
point of the domain. If the domain is terminated at an inner boundary 
p = p (corresponding to r = b), then the coefficients B m may be different 
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from zero, and their values would have to be determined by resorting to an 
appropriate numerical technique, such as approximation in the mean, together 
with an iterative scheme . 


THE SOLUTION I FOR THE PARACHUTE TYPE STRUCTURE 


For the parachute type structure the stress function I must satisfy 
equation ( 36 ) subject to the condition ( 37 ) • The solution for this case is 
similar to that for the single scallop and can be written as 


¥(p,x) = — p 2 ^ a^x 11 

n=o, 2 , 4 , . . . 


+ HP* 


Y Pn* 11 - log(p 41 + x 2 ) 


n= 0 , 2 , 4 , . . 


Y fmW + B mP Am S m W 


(A42) 


m 1 , 2 , 3 , • • • 


The coefficients a n , P n and the functions f m , g m can be found using the 
approach outlined earlier in this section, where p a and are expressed 
as follows : 


3(3 , 2 t A n 

- m a ^) x o 


Ma = 


- 1 - L (?K - Vt 

_ n=o,2,4^. > A n + ^ 

L ( 3 < - ¥^1 ^ 

n=o, 2 , 4 , ...V n+1 n+1 J 


(a4 3 ) 


. 2 . 5 ,2 . 1 v4 . V /Cpp' _ 4 + 5(3 p» 4p n- 2 Vn 
3 9^ 5 x o n=0 b ^rnr^ wtt) 


“e ■ - 


x r 


n= 0,2,4, 


ri 4+58 „ 4p" 

5pa n - - - - ? - p - p" a 

n n+1 11 n + 1/ 0 




(A44) 


The value of the constant k in (a 42) depends on the tension in the ropes in 
the region p s p Q . 

For the parachute type structure it can be verified that the function 
u^x) must satisfy the boundary condition 

(A^ - l)x 0 u m (xo) + (1 + x|)u^(x 0 ) = 0 (A45) 


where 
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m 


^ ^ r x 0 UmUm dx 

•To (TT^- 


~ X + °^2 %( X 0 )%( S 0^ 



(1 + X 2 ) 2 . 


Um% dx 


Xq 

1 + x o 


u m( x o) u m( x o^ 


(A46) 


It follows from equation (A 46) that A^ must always be real. The positive 
definite character of A 2 is difficult to establish in this case, but since 
the solutions corresponding to negative A 2 have a rapidly oscillating 
singularity as p becomes small, the undetermined coefficients associated 
with these solutions may be chosen appropriately as zero. The eigenvalues of 
interest, then, are related to the positive values of Aj| and are easily 
obtained by the interpolation procedure described earlier in this appendix. 
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APPENDIX B 


MEAN SQUARE APPROXIMATION FOR COEFFICIENTS A m 


The method of mean square approximation, used in determining the coeffi- 
cients A m of equations ( 28 ) and (4l), will he described briefly in this 
appendix. Consider the case of the single scallop subject to the condition 


1 8 ^ J3_ 

_P 8p p 2 ^x 2 



at p = 1 


(Bl) 


Using 1, as defined in (28), the above condition can be expressed as 


I 

m=i,a,3 , 



. n=o, 2,4 


+ + 2)(n + 



n=0,2,4 


[ 2 %i + P(n + 2 )(n + lja^+^jx 13. 


+ p ^ l>P n + P(n + 2)(n + l)(3 n+ 2 ]x ri - + x 2 ^ = 0 


n .~0 ^2 ^4 ; • ■ * 
or, more concisely, as 


2 A mSm( x ) " yW = 0 


in 1 , 2 , 3 , • • 


where 


Sm(x) = zL 


a x 11 
m,n 


n=o,2,4 


? * * • 


y(x) = ^ u n x n 

n=o, 2 , 4 , . . . 

The coefficients a m ^ n , in equation (b 4) are defined by 


(B2) 


(B3) 


0*) 
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where 


^l^n ^eAujII + P( n + 2)(n + l)C m ^ n+2 

Mn = - Jjj [ 2a ti + P(n + 2)(n + 1)%^] 

- M-t^Pn + P( n + 2)(n + l)Pn +2 - & n ]^ 

& ft = - , 8 =1 
0 3 2 

6 n = 0 for n > 2 , 


(B5) 


(b6) 


The summation in (B3) is an infinite series, which includes correspond- 
ingly an infinite number of the undetermined coefficients A m . For practical 
purposes this series is truncated after a suitable number of terms, k, and 
the coefficients A m are determined to minimize the square of the difference 

between the sum E A m g m (x) and the function y(x) in the interval 

k 

0 s x s: Xq . This is the basic principle underlying the technique of mean 
square approximation (ref. 11 ) . Its application to scalloped paraboloids can 
be described as follows: 


Define 


and 


R(x) = 


I 


'm 


(x) - y(x) 


. . .k 

> x of 


Q - / [R(x)r dx = 

J o J o 


r i 


A m g m (x) - y(x) 


dx 


un=i, . . .k 

For Q to be minimum the following conditions must be satisfied: 

d 


dAr 

Hence , the condition 


(Q) = 0 for r = 1, 2, . . . , k 


(B7) 


(B8) 


(B9) 


/ [2 A m g m (x)g r (x)]dx = / y(x)g r (x)dx r = 1, 2, . . . k 

J o 

or, interchanging the order of integration and summation. 


(BIO) 


Z p x o p x o 

A m J g m ( x )s r (x)dx = J y(x)g r (x)dx 


r = 1, 2, 


m=i,2. 


k 

(Bll) 


3 A 




y(x)g r (x)dx 


for r,m = 1, 2, . . . k 


x 0 S ' m+1 (Bl4) 

n + s + 1 ° 


S — O ,2 ,4 , • • « 11— O ,2 ,4 , • * • 


The unknown coefficients A m are determined by solving the simultaneous 
equations (B12) . 


The results presented for the two examples have been obtained by the 
technique just described. When the structure has a boundary at p = p g 
(corresponding to r = b) and the magnitude of the lowest eigenvalue, [A-J , 
is sufficiently high, the homogeneous solutions are such that any disturbance 
produced at the edge p = 1 decays before reaching the edge p = p Q . The 
method given here, then, can be used separately for correcting the value at 
each of the two boundaries, p = p Q and p = 1. In other cases it may be 
necessary to combine the present approach with an iterative scheme similar to 
the one described in reference 8. 
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Figure 4.- The single scallop; membrane force n r along the edge; p 

0 S x s x^. 





Figure 6.- The single scallop; membrane forces n r , n# along the central 
meridian; x = 0; 0 £ p s; 1; n r g = 0. 
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Figure 12.- Parachute type structure; membrane forces n r , n q , n r g along the 

edge; x = p = 0.1 £ p :£ 1. 
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Case H : Stresses at 9=0 
h = 0.4 a 
N = 12 
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Figure 15 •- Parachute type structure; membrane forces n r , ng along the 
central meridian due to the solution for selected values of |3; 

p Q £ p = 1; n r g =0. 
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